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ABSTRACT 

Investigation   of  oceanographic  data   acquired   in   the  Columbia    River 
estuary,   which   lies  between  the  states  of  Washington  and  Oregon, 
revealed   significant  variation   in   sound   velocity  with    respect  to  time  and 
location.      National   Ocean   Service   (NOS)    standards  for  echo-sounder 
acquired   soundings    require  the   knowledge  of   sound   velocity  to  within 
plus  or  minus  4  meters   per  second.      To  meet  this    requirement   in   the 
Columbia    River  estuary,    methods   were  devised  to  allow   prediction   of 
sound   velocity   based  on    location   and  the   height  of  tide. 

Confidence   intervals   associated   with    sound   velocity   predictions 
provided   by   regression   analyses   showed   a   substantial    improvement   in 
accuracy  when   compared   to   using   a   single  average  velocity   for  each 
location.      Nevertheless,    there  were  still    some  cases   where  the   NOS 
requirement  could   not  be  met. 
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I.  INTRODUCTION 

A.   NOS  HYDROGRAPHY 

"Hydrography   is  that  branch   of  physical   oceanography  dealing   with 
the  measurement  and   definition  of  the  configuration  of  the  bottoms   and 
adjacent   land   areas  of  oceans,    lakes,    rivers,    harbors,    and  other  water 
forms  on   earth"    [Ref.    1:    p.    1-3].      The   National   Ocean    Service    (NOS) 
is  the  federal   agency    responsible  for  producing   accurate   nautical   charts 
of   nearshore  areas  of  the  United   States   and   its  territorial   possesions. 
The  principal   objective  of  most   hydrographic  surveys   conducted   by   NOS 
is   to  obtain   basic  data  for  the  compilation   of   nautical   charts   with 
emphasis   on   the  features   that  may  affect  safe   navigation    [Ref.    1:    p. 
1-3].      For  this  thesis   data   analysis   is    limited  to  the  study  of  sound 
velocity  variations   and  the  correction   of  echo   soundings   for   NOS 
hydrographic   surveys.      The  type  of  data   used  for  this   investigation    is 
often   the  basis   for  other  physical   studies   such   as   the  dynamics  of 
oceanic  mixing   and   salt  transport.      No  attempt  will   be  made  to  analyze 
or  explain    such    processes. 

1 .  Sound  Velocity  and  Sound  Speed 

Manuals   and   papers   about   hydrographic   surveying   often    use 
the  term   "sound   velocity"   to  designate  what  a   physicist  would   define   as 
"sound   speed."      The  term   "sound   velocity"   will    be   used    in   this   thesis 
as   many  of  the   references   also   use  the  term. 

2.  Echo  Sounding 

Echo   sounders   are  the   basic  tool   of  the   NOS    hydrographer  for 
determination   of  depth.      Echo   sounders   measure  the  time   required   for  a 
sound  wave  to  travel   from   its   point  of  origin   to  the  bottom  and    return. 
They  also  convert  this   measured   time  to  distance  or  depth    [Ref.    1:    p. 
AG-1].      The   basic  formula   for  echo   sounder  depth    determination    is: 


Z  =     VT/2  (l.D 

where  Z  is  the  depth,   V  is  the  calibrated  velocity  of  sound  for  the  echo 
sounder,    and  T  is  the  time  between  transmission  and   reception  of  the 
sound   signal.      For  depths  to  be  accurate  precise  measurements  of  time 
and  a   knowledge  of  the  mean,   or  harmonic,    sound  velocity   (Chapter   II) 
are  mandatory. 

B.  THE  VELOCITY  OF  SOUND   IN   SEA  WATER 

Wilson    [Ref.    2:    p.    641]    developed   empirical   equations   based  on 
laboratory  measurements  of  sound   velocity   under  tightly  controlled 
conditions.      The  basic  equation   is: 

V  =   1449.14  +  VT   +  Vp  +  Vs   +  VSTp  (1.2) 

where  V  is   the  velocity  of  sound   in   meters   per   second,    V-p   is   a 
correction   for  temperature,    Vp   is   a   correction   for  pressure,    V^    is   a 
correction   for  salinity,    and  V^jp   is   a   simultaneous   correction   for  the 
combined   effects  of  temperature,    salinity,    and   depth    [Ref.    1:    p.    4-75]. 
Wilson's   equations    (Chapter    II),    which   are  currently    used   by   NOS, 
were  chosen   to  analyze  the  data   for  this   thesis. 

C.  NOS  ACCURACY   REQUIREMENTS 

The  sound   velocity  at  which    NOS   echo  sounders   are  calibrated    is 
800  fathoms   per  second.      To  obtain   accurate  depths   for   nautical   charts 
the   hydrographer  must  correct  all    soundings   where  the  velocity   of 
sound   deviates   from   the  calibrated   value.      There  are  several   methods 
available  for  determining   corrections   for  soundings,    and   they  will    be 
discussed   in   the  following   section.      The  NOS   accuracy    requirements 
[Ref.    1:    p.    4-70]    are: 


For   use   in   correcting   echo  soundings,    the  velocity  of   sound   must 
be    known   with    sufficient  accuracy   to  ensure  that   no   sounding   will 
be   in   error  by   as   much   as   0.25   percent  of  the  depth   from  this 
cause  alone.      Therefore,    the  mean   velocity   must   be   known   to 
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within   plus  or  minus  4  meters   per  second.      To  satisfy  these 
accuracy   requirements,    one  must   know  the  mean   temperature  of  the 
water  to  an   accuracy  of  plus  or  minus   1    degree  Celsius   and  the 
salinity  to  within   plus  or  minus   1    part  per  thousand.      A   sufficient 
number  of  salinity  and  temperature  observations   must  be  made  so 
that  the  velocity  can   be  determined  within  the  specified  accuracy 
over  the  entire  area   sounded. 


D.   PROCEDURES  FOR  CORRECTING  SOUNDINGS 

Echo  soundings   can   be  corrected   by  either  a   direct  comparison   or 
an   oceanographic  procedure.      Direct  comparisons   involve  comparing   echo 
soundings  to   known   depths;    this   does   not   require  a    knowledge  of  the 
sound   velocity.      Procedures   which    involve  the  calculation   of  sound 
velocity  from  measurements   of  water  characteristics   are  classified   as 
"oceanographic  and    limnologic  determinations"    [Ref.    1:    p.    4-74].      For 
this   thesis  the  direct  measurement  of  sound   velocity   by  a   sound 
velocimeter  is   also  classified   as   an   oceanographic   procedure. 

1  .        Direct  Comparison  Systems     * 

Two   procedures    used   to  determine  depth    corrections   by   the 
direct  comparison   method   are  the  bar  check   and   the  vertical   cast.      In 
these  the  effects  of   sound   velocity   variation   and    instrument   error  are 
included   in   the   results    [Ref.    1:    p.    4-71]. 

a.        Bar  Check 

The  bar  check   procedure   involves    lowering   a    bar  over 
the   side  of  a    survey   vessel    and   positioning   the   bar   under  the 
transducer  at  a    known    depth    indicated   by   markings   on   the   lines 
supporting   the   bar.      Comparisons   between    the  depth    as    indicated   by 
the  echo   sounder  and   the   known    depth   are  made  at   5-   or  6-foot 
increments   as   the   bar   is    lowered   and    raised   though   the  water  column. 
Due  to  the   handling   of  the   bar,    this   procedure   is    usually   limited   to 
smaller  vessels   and   the  quality  of  the   results    is   a   function   of   sea 
conditions.      Observations   during    rough    seas   are   subject  to   unacceptable 
magnitudes   of  error    [Ref.    1:    p.    4-71]. 
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b.       Vertical   Cast 

The  vertical  cast  requires  simultaneous  sounding  of  the 
bottom  using  an  echo  sounder  and  a   lead   line.      This  procedure,    which 
can   be  accomplished  by  practically  all  survey  vessels,    is   not  as 
accurate  as  the  bar  check,    but  vertical  casts  can   be  used  to  determine 
or  verify   instrument  errors.     Where  sounding  sheets  junction,    there  are 
often   discrepancies   and   soundings  fail   to  agree.      Vertical   casts   can 
sometimes   help   resolve  such   problems    [Ref.    1:   p.   4-72]. 

2.        Oceanographic  Sound  Velocity  Determinations 

The  more  traditional   oceanographic  procedures   involve  sound 
velocity  determination   by  measurement  of  the  variables   which   enter 
Equation    1.2,    namely  temperature  and   salinity  as   functions  of  depth. 
The   Nansen   cast  and   electronic   sensors   such   as   the  CTD    (conductivity, 
temperature,    and   depth)   or  the  STD    (salinity,    temperature,    and   depth) 
are  examples   of  equipment   used   for  these  procedures.      An   alternative   is 
the  direct  measurement  of  sdund   velocity   using   a   sound   velocimeter. 
The  electronic   sensors    record   data   almost  continuously  with   depth. 
[Ref.    1:    p.   4-74]. 

To  obtain   velocity  corrections   with    Nansen,    STD,    or  CTD 
data,    Umbach    [Ref.    1:    p.    4-75]    recommends   using   the   "summation   of 
layers"   method,    which   extends   from  the  concept  of  mean   or   harmonic 
sound   velocity. 

Most  of  the  electronic   sensors   which   are  currently   in    use   yield 
precise  measurements    in    excess   of   NOS   accuracy    requirements   for   sound 
velocity  corrections    if  properly   calibrated   and   maintained    [Ref.    1:    p. 
4-74].      According   to   NOS   standards   all   electronic   sensors   must  be 
periodically   checked   for  accuracy   against   Nansen   cast  or  equivalent 
observations. 

E.       SOUND   VELOCITY   VARIATION 

The  effect  of  changes    in   temperature,    salinity,    and   pressure  on 
sound   velocity  depends   upon   the   initial   conditions  of  these  variables. 
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Wilson's   equations  are  used  to  illustrate  in   Table   I   the  effect  of 
changing   variables   by  one   unit. 


TABLE   1 

Variation  of  Sound 

Velocity 

PRES 
db 

TEMP 
°C 

SAL 
g/kg 

V 
m/s 

delta   V 
m/s 

1.00 
1.00 

15.00 
14.00 

30.00 
30.00 

1502.90 
1499.30 

3.60 

1.00 
1.00 

5.00 
4.00 

30.00 
30.00 

1464.27 
1460.03 

4.25 

1.00 
1.00 

10.00 
10.00 

35.00 
34.00 

1490.76 
1489.48 

1.28 

1.00 
1.00 

10.00 
10.00 

1.00 
0.00 

1449.01 
1447.84 

1.17 

60.00 
59.00 

10.00 
10.00 

30.00 
30.00 

1485.33 
1485.32 

0.02 

1.00 

o.oo. 

10.00 
10.00 

30.00 
30.00 

1484.38 
1484.36 

0.02 

Velocity    (V)    was    c 

omputed  usi 

rig  Wilson' 

s    equations. 

Corrections   for  .sound   velocity   variations    in   estuaries   can    present 
the   hydrographer  with   a   difficult  task.      According   to   Cameron    and 
Pritchard    [Ref.    3:    p.    306],    an    estuary   is   a   semi-enclosed   coastal    body 
of  water   having   a   free  connection   with   the  open    sea   and   within   which 
the   sea   water   is   measurably   diluted   with   fresh   water  deriving   from   land 
drainage.      Estuaries   are  extremely   dynamic   systems.      Topography, 
river  flow,    and   tidal   action   are   important  factors   that   influence  the   rate 
and   extent  of  mixing   salt  and  fresh   water   [Ref.    4:    p.    4].      For  the 
hydrographer,    estuaries   are  areas   where   sound   velocity   varies    not  only 
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with   location   but  also  in  time  due  to  the  tidal  cycle  and  variations   in 
river  flow.      Although   no  specific  procedures  are  given,    Umbach 
[Ref.    1:   p.   4-75]    notes  the  possible  complexity  of  such  areas: 


Under  special  circumstances  such  as    ...    in  areas  of  extensive 
estuarine  discharge  of  fresh   water   .    .    .    series  of   regional   vejocity 
curves   are   required.      Regional   curves   must  be  carefully   studied   to 
determine  how  they  can   be  best  grouped  by  area  or  time  or  to 
permit  drawing  of  average  curves. 


F.       OBJECTIVES 

One  objective  of  this  thesis  was  to  investigate  the  magnitude  of 
sound  velocity  variation  associated  with  significant  changes  in  water 
characteristics   within   an   estuary. 

A  second  objective  was  to  investigate  the  feasibility  of  determining 
velocity  corrections  according  to  NOS  standards  in  areas  where  extreme 
changes   in   water  characteristics  occur  with    respect  to   space  and   time. 

A  third  objective  was   to  write  a    FORTRAN    program  to  enable 
further  testing   by  the  field   hydrographer  equipped   with   a   personal 
computer.      The   use  of  commercial    subroutines    has   been   avoided. 
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II.  METHODOLOGY  AND  PROCEDURE 

A.  AREA  OF  STUDY 

The  Columbia    River  estuary  which   discharges   into  the   Pacific 
Ocean   at  latitude  46°15'N,    longitude   124°00'W  has   variations   in   water 
characteristics   which   make   it  an   excellent  area  for   investigation. 
According  to   Dyer   [Ref.    4:    p.    96] ,    the  Columbia    River  estuary   is 
characterized   by  an    unusual   combination   of   large   river  discharge  and 
strong  tidal   currents.      Recent   investigations    have   revealed   differences 
in   salinity  of  up  to  30  grams   per   kilogram  at   individual   stations   during 
a   single  tidal   cycle    [Ref.    5:    p.    104].      Due  to  this    large  difference  and 
the  availability  of  data  from  a    recently  completed   survey,    this    region 
was   chosen   for  investigation. 

B.  DATA   SOURCES 

To  analyze  the  problem  completely,    water  characteristics 
(temperature  and   salinity   as   functions   of  depth)    must   be   known   at 
stations    spaced   throughout  the   region,    and   each    station   occupied   during 
a   complete   lunar  tidal   cycle.      Variations   caused   by   seasonal   changes    in 
fresh-water   runoff  should   be  examined   thoroughly. 

1  .        Columbia  River   Estuary  Data  Development  Program 

The   Columbia    River   Estuary   Data    Development   Program 
(CREDDP)    was   a   5-year,    federally   funded   study   which    included 
approximately   2   years   of  data   acquisition   for   investigation   of   physical 
and   biological   characteristics.      CREDDP  was   designed   to  meet  the   needs 
of  two  major  groups.      One  group   consisted  of   local,    state,    and   federal 
government  agencies   involved   in   planning   and   permitting   activities. 
The   second   group   included    research    scientists   and   educational 
institutions    [Ref.    6:    p.    v] . 
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One  accomplishment  of  CREDDP  was  the  acquisition  of  current 
velocity,    conductivity,   temperature,   and  depth    (VCTD)   data  for  a 
better  understanding  of  physical   processes  such   as  currents,    sediment 
transport,   and  salinity  intrusion.      VCTD  profiling  was  carried  out  for 
13  consecutive  days  from  17  stations  for  periods  of  from  12  to  48  hours 
per  station.      Between   October  16  and  27,    1980,    approximately   1600 
VCTD  casts  were  made.      Of  these,   458  casts  were  made  available  to 
users    [Ref.   6:    p.    L-8] . 

Data  were  supplied  to  NPS   by  the   U.S.    Army   Corps   of 
Engineers    (COE)    at   Portland,    Oregon,    on   two  9-track  magnetic  tapes 
[Ref.   6:    pp.    1-5]. 

For  an   analysis   of  sound   velocity   variation,    VCTD   casts   from 
11    of  the  17   stations   were   investigated.      The  stations,    along   with   their 
positions   and   times   are  listed   in   Table   II.      The  plot  of  the  geographical 
positions   of  the  stations   is   shown    in    Figure  2.1.      These   stations 
covered  the  area  where  the  maximum   sound   velocity   variation   would   be 
expected   as   deduced   from   salinity   intrusion   diagrams  of  the  estuary 
[Ref.   5:    p.    104]. 

2.  NOS  Bar  Check  Data 

NOS   bar  check  data  were   used   for  a   sound   velocity   variation 
analysis.      The  descriptive   report  for  survey   H-8422   listed   the  day   each 
bar  check  was  taken   and  the  correction   for  each   5-foot   increment  of 
depth.      Also  given   was   a   general   term  for  the  stage  of  the  tide   such   as 
ebb,    flood,    high   water,    or   low  water.      Data   for  46   bar  checks   were 
acquired   between   May  9   and  July   16,    1958    [Ref.    7:    p.    1]. 

3.  NOS  Predicted   Tide  Data 


NOS   predicted   tide  data   were  obtained   for  correlation   and 
regression   with   the  CREDDP   and  the   H-8422  bar  check  data.      Predicted 
tides   were   used   because  field   hydrographers    rarely   have  the  time  or 
the  proper  tools   for  a   detailed   analysis  of   real   tide  data. 
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TABLE 

II 

Columbia  River 

Estuary  Data  Development  Program 

STATION  (STN)  INFORMATION 

STN 

POSITION 

ACQUISITION 
START 

TIMES  (PDT) 
END 

2N 

46015'45"N 
124°03"31"W 

10/18/80 

0930 

10/19/80 

0230 

2S 

46015'13"N 
124°03*20"W 

10/18/80 

1230 

10/19/80 

0930 

4NA 

46°14'03"N 
123°56'09"W 

10/20/80 

2115 

10/21/80 

1745 

4NB 

46°14'47',N 
123°55*53"W 

10/20/80 

2105 

10/21/80 

1730 

4SA 

46012'42'"N 
123°57'09,,W 

10/26/80 

0007 

10/27/80 

0832 

4SB 

46°12'51,,N 
123°56'55"W 

10/25/80 

2350 

10/27/80 

0820 

5NA 

46°13'55"N 
123°53' 17"W 

10/20/80 

1733 

10/22/80 

0930 

5NB 

46°14' 18"N 
123°53' 14"W 

10/20/80 

2035 

10/22/80 

0913 

5NC 

46°14'07"N 
123°53' 14"W 

10/21/80 

0935 

10/22/80 

1003 

5SA 

46°11*55"N 
123°48'38"W 

10/24/80 

0830 

10/25/80 

233  1 

5SB 

46°12'05"N 
123°48'54,,W 

10/24/80 

1126 

10/25/80 

1130 
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Figure  2.1        CREDDP   Stations. 

C.       CREDDP   DATA   ANALYSIS 

1  .        COE   VCTD   Data   Tapes 

COE  VCTD   data   tapes   contained   pertinent   information    such    as 
header  data,    time,    temperature,    depth,    conductivity,    and    salinity. 
Mercier    [Ref.    6:    p.    L-8]    shows   a    breakdown   of   the   data   components   for 
one   record   of   data.      The   tapes   contained   more   data    records   than    were 
necessary   for   this    investigation.      The  VCTD   sampling   device   recorded 
data   every   0.04   seconds,    thereby   making   the  average   number  of 


18 


records  for  a  single  cast  approximately  1500.      Having  270  casts  to 
analyze,    a  method  was  devised  to  reduce  the  data  to  a  more  efficient 
size.      A  computer  program  was  formulated  to  input  the  data  for  each 
cast  and  then  output  the  first  record  and  the  records  which 
corresponded  to  the  nearest  whole  meter  of  depth.      For  example,    a 
10-meter  cast  would  theoretically   have  11    records  of  data;    the  first 
record  would    represent  the  surface  layer,    and  the   remaining    10   records 
would   be  for  1-meter  increments.      In   accordance  with   the   NOS   policy  of 
using   bar  check  data   measured   at  5-   or  6-foot   increments,    data  were 
eliminated   except  at  1 -meter   increments.      A   sample  of  a   profile  after 
editing   is   shown    in   Table   III.      Standard   depths   which    represent  depths 
of  1,3,    5,    7,    9,    11,    and   13  meters,    were  adopted  to  further  simplify 
data   processing   and   description. 

Several   drawbacks   were  encountered   using  the  CREDDP   data. 
Because  of  the  shallow  depths  of  CREDDP   stations,    many  of  the  casts 
were   inadequate  for  correcting   soundings   in   the  deeper  areas   of  the 
estuary.      Some  casts   did    not   include   surface   layer  data    (0  to  2 
meters) . 

2.        Sound  Velocity  Computations 

a.        Wilson's    Equations 

Sound   velocity  was   computed  for  each    record   by  Wilson's 
equations   which    use  temperature,    pressure,    and   salinity   for  the 
computation.       In   the  following   equations   V   is   sound   velocity    in   meters 
per   second;    T   is   temperature   in    degrees   Celsius,    where   -4   <   T   <  30 
°C;    P   is   pressure   in    kilograms   per   square  centimeter,    where   1    <    P   < 
1000    kg/cm    '    and   S    is    salinity   in   grams   per   kilogram,    where  0  <   S   <  37 
g/kg. 

V  =  1449.14  +  VT     +  Vp     +  Vs     +  VSTp  (2.1) 

VT     =  4.5721T   -  4.4532x10~2T2   -  (2.2) 

2.6045x10"3T3   +   7.9851x10"6T4 
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VP 


=  1.60272x10_lP  +  1.0268x10"5P2  ♦  (2.3) 

3.5216x10"9P3  -  3.3603x10"12P4 


Vs     =  1.39799(S-35)   ♦  1 .69202x10"3(S-35)2  (2.4) 

VSTp       =   (S-35)(-1.1244x10"2T  +  7.7711x10~7T2  ♦  (2.5) 

7.7016x10~5P  -   1.2943x10"7P2  ♦ 

3.1580x10"8PT  +   1.5790x10~9PT2)    ♦ 

P(-1.8607x10'4T  +  7.4812x10~6T2  * 

4. 5283x1 0"8T3)   +   P2(-2.5294xl0"7T  + 

l.8563x10"9T2)    +   P3(-l.9646x10"10T) 
The  standard   deviation   for  the  velocity   is  0.30  meters   per 
second   [Ref.   8:   p.    1357]. 


TABLE 

III 

4SA 

--Cast  6 

--Sound 

Velocity  Variables 

COLUMBIA   RIVER 

STATION   = 

4SA 

CAST    =    031B08 

LATITUDE 

46°12' 

42  "N 

LONGITUDE    123°57'09"W 

DATE    OF    CAST    1980 

OCT    26 

TIME    OF    CAST       5: 16:     0 

CTD 

TIME 

PRES 

DEPTH 

TEMP      COND 

SAL 

s 

db 

m 

°C      ratio 

g/kg 

0.04 

-0.04 

0.  10 

12.26    0.4531 

15.71 

5.48 

0.85 

1.00 

12.04   0.5252 

18.57 

12.04 

1.86 

2.01 

11.79    0.7073 

25.91 

18.48 

2.85 

3.00 

11.78    0.7571 

27.95 

24.60 

3.86 

4.01 

11.78    0.7932 

29.42 

31.52 

4.86 

5.00 

11.79    0.8184 

30.46 

37.68 

5.90 

6.04 

11.79    0.8164 

30.37 

43.80 

6.86 

7.00 

11.80    0.8232 

30.64 

50.48 

7.88 

8.01 

11.81    0.8291 

30.88 

56.80 

8.89 

9.01 

11.81    0.8302 

30.92 

63.32 

9.88 

10.00 

11.81    0.8330 

31.04 

69.92 

10.91 

11.03 

11.82    0.8338 

31.07 

76.12 

11.90 

12.01 

11.82    0.8356 

31.  14 

82.84 

12.90 

13.01 

11.82    0.8346 

31.10 
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b.       Mean   Sound  Velocity 

Because  echo  sounders   are  time  measuring   devices, 
distance  or  depth   is  one-half  the  round  trip  travel  time  of  a   sound 
pulse  multiplied   by  a   suitable  mean   velocity.      The  mean   sound   velocity 
(MSV)    is  given   by 

MSV  =  Z/T  (2.6) 

where  Z   is  the  distance  traveled   in   time  T.      To  correct  for  sound 
velocity  variations,    the  water  column   is   divided   into  a   series   of   layers, 
each   with   thickness  AZj   and  an   associated   velocity  Vj .      Therefore,    the 
time   interval   for  the  sound  wave  to  pass   through   the   i        layer   is: 

ATj     =  AZ/Vj  (2.7) 

Summing   all    layers    results   in: 

T     =    2     AT:        =2     AZ:/V:  (2.8) 

1-1     '      1=1     '    ' 

or: 

T   =  /gdZ/V  (2.9) 

Substituting   into   Equation   2.6   results   in: 

MSV  =  Z/(/*dZ/V)  (2.10) 

Equation   2.10  was   evaluated    using   the  trapezoidal   method    [Ref.    9:    p. 
86].      NOS    uses   mean    sound   velocities    in   the   "summation   of   layers" 
method   described   by   Umbach    [Ref.    1:    p.    4-74]    and  Wallace    [Ref.    10: 
p.    RK530.1]. 

Two   sound   velocity   listings   computed   using   the  data    in 
Table    III    are   shown    in   Table    IV.      The   sound   velocity   column   was 
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computed   using   Equation  2.1,    and  the  mean  velocity  column  was 
computed   using   Equation  2.10. 


TABLE   IV 

• 

4SA--CAST 

6--Sound  V 

elocity  Data 

DATE: 

26    OCTOBER    1980 

TIME: 

0520 
SOUND 

(PDT) 
MEAN 

DEPTH 

VELOCITY 

VELOCITY 

m 

m/s 

m/s 

0.10 

1475.5 

1475.5 

1.00 

1478.1 

1476.8 

2.01 

1486.1 

1479.6 

•     3.00 

1488.6 

1482.2 

4.01 

1490.5 

1484.1 

5.00 

1491.8 

1485.5 

6.04 

1491.7 

1486.6 

7.00 

1492.1 

1487.3 

8.01 

1492.5 

1488.0 

9.01 

1492.5 

1488.4 

10.00 

1492.7 

1488.9 

• 

11.03 

1492.8 

1489.3 

12.01 

1492.9 

1489.5 

13.01 

1492.9 

1489.8 

14.01 

1492.9 

1490.0 

To   simplify  the  analysis   of  the   variation   of  mean    sound 
velocity  as   a  function   of  time,    a   computer   program  was   written   to 
arrange  the  data   in    individual   arrays   for  each   depth   at  each    station. 

3.        Mean  Sound  Velocity  and   Tide  Height 

Inspection   of  the  mean   sound   velocity   arrays   suggested   that 
velocities   were  maximum   at   high   tide  and   minimum  at   low  tide.      To  test 
this,    a   correlation   and    regression   analysis   between   the  CREDDP  derived 
mean    sound   velocities   and   predicted   tide   heights   was   performed. 
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a.        Predicted  Tide  Heights 

To  perform  the  correlation  and   regression  analysis,    a 
method   had  to  be  developed  to  calculate  the  height  of  the  tide  at  any 
desired  time.      One  method  would   be  to  obtain   the  harmonic  constants 
for  prediction   of  the  Columbia    River  estuary  tides  from  the  NOS   Tides 
and  Water   Levels    Branch.      These  constants  would  then   be  used    in   the 
appropriate  equations    [Ref.    11:    p.    123].      From  a   computer  programming 
standpoint  this  method  would   be  the  simplest.      A   second   method   would 
be  to  formulate  a   computer  program  based  on   the  graphical   method 
(Figure  2.2)    described   in   the   NOS   Hydrographic  Manual    [Ref.    1:    p. 
4-68].      This   method  was   chosen   because   normally  the  field 
hydrographer  is   not   supplied   with   the   harmonic  constants   for  project 
areas.      However,    the  tide  tables   are  an   essential   tool   of  the 
hydrographer. 

The  present   NOS   HYDROPLOT  computer  program  for  tide 
prediction   lacks   documentation   and   is   not  designed   to   solve  this 
particular  problem    [Ref.    10:    p.    500.1].      Therefore,    a   program  was 
written. 

Points   A   and    E    (Figure  2.2)    represent  the   predicted    low 
and    high   tide,    respectively.      The   straight   line   between   these   points 
was   divided    into  four  equal    parts   by   points    B,    C,    and    D.      Point   B'   was 
plotted   directly   below   B   at  one-tenth   the  total    range,    and   point   D'   was 
plotted   directly   above   D   using  the  same  distance   [Ref.    1:    p.    4-68], 
forming   three   regions   for  the   solution   of   predicted    heights   between 
points   A   and    E.      For  the   region    between    points   A   and    B',    the  curve 
was   assumed   to   be   parabolic,    and    in   order  to  calculate  the  tide   height 
for   a   designated   time   in   this    region,    the   equation 

Hj     =  aTj2   +   bTj    +   c  (2.11) 

was    used,    where   H    is   the  tide   height  and   T   is   the  designated   time.      To 
calculate  the  coefficients   a,    b,    and   c;    the  equations 

HA     =  aTA2   +   bTA   +   c  (2.12) 
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Figure  2.2        Predicted   Tide  Graph 
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HB'     =  aTB'2  +  bTB'   *  c  (2-13^ 

0     =  2aTA  +  b  (2.14) 

were  used.      For  times  that  lie  between   points   D'   and   E,    the  same 
method  was   used.      Between   points    B'   and   D'   the  change   in    height  was 
considered  to  be   linear  with    respect  to  time  and   can   be  computed   by 
the  equation 

Hj     =  HB.   ♦  m(TrTB0  (2.15) 

where  m   is   the   slope  defined   by 

m  =    (HD--HB-)/(TD--TB0  (2.16) 

Using  the  program   based  on   the  above  equations,    tide   heights   were 
calculated   for  the  times  of  all  VCTD  casts   used   in   the  analysis.      "Data 
Array"   means   a   set  of  sound   velocity   and   predicted   tide   height  data 
that   shows   the   variation   of  both   variables   with   time  for  a   particular 
station   and   depth    (Table  V) . 

b.  Data    Plots 

Computer  generated   plots   of   predicted  tide   height  and 
observed   mean    sound   velocity    (Figure  2.3)    were  made  for  all    stations 
using   mean    velocities   based   on    standard   depths. 

c.  Correlation   Analysis 

A   correlation   analysis    routine  for  mean    sound   velocity 
and   tide   height  was    run   for  all   data   arrays.       In   addition   to  computing 
the  correlation   coefficient,    the   routine  also  computed   the  averages   and 
the   standard   deviations  of  the  mean    sound   velocity   and   the   predicted 
tide   height  for   each   data   array.      Averages   were  computed   by  the 
equation 

X  =     SxVn  (2.17) 

1  =  1    ' 
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TABLE 

V 

4SB--Data 

Array-- 

•Depth  =  7  m 

TIDE 

MEAN 

TIME 

HEIGHT 

VELOCITY 

hr 

ft 

m/s 

0.12 

1.2 

1469.6 

1.50 

4.1 

1476.6 

2.63 

6.5 

1486.5 

3.60 

7.6 

1490.8 

4.62 

7.6 

1490.8 

5.27 

7.1 

1487.3 

6.32 

5.6 

1483.1 

9.23 

1.7 

1472.1 

11.55 

3.0 

1471.0 

12.58 

4.9 

1472.2 

13.17 

6.0 

1476.9 

17.53 

7.7 

1489.6 

21.52 

-0.8 

1471.4 

22.43 

-1.5 

1468.6 

23.35 

-1.3 

1467.7 

24.28 

-0.2 

1467.1 

25.33 

1.9 

1470.0 

26.42 

4.1 

1473.6 

27.57 

6.2 

1480.0 

28.48 

7.2 

1485.9 

29.48 

7.2 

1491.9 

30.50 

6.4 

1488.6 

31.53 

5.1 

1485.4 

32.53 

3.8 

1477.3 

J.L. 

where   X:    is   the   i        observation   and   n    is   the   number  of  observations 
[Ref.    12:    p.    124].      Standard   deviations    [Ref.    12:    p.    128]    were 
computed    using   the  formula 


s   =    rk    (xrx)2/(n-1)  (2.18) 

v  i  =  l 

Computation   of  the  simple  correlation   coefficient    [Ref.    12:    p.    258]    of 
the  mean    sound   velocity   and   predicted   tide   height  arrays   was   done  with 
equation 

r=  [n  n4  XjY,    -    (JfrH-l^  Yj)] (2.19) 

I       n      z  n z n ~ n ^ — 

/[n  2    Xj2   -    (.2    Xj)2][n.2    Y;2-    (   2    Yj2)] 

V  1=1     ■  i=l     ■  i=l      '  i=i      ' 
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Figure  2.3       4SB--Data   Array   Plot--Depth   =   7  m. 
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A  sample  of  the  results  of  a  correlation  analysis   is   listed   in   Table  VI 


TABLE  VI 
4SB--Statistical  Data—Depth  =  7  m 


AVERAGE    TIDE   HEIGHT    4. 4    ft 

STANDARD   DEVIATION   OF   TIDE    HEIGHT    3 . 0    f t 

AVERAGE    MEAN   VELOCITY    1480 . 8    m/s 

STANDARD   DEVIATION   OF   MEAN   VELOCITY    8.0    m/s 

CORRELATION   COEFFICIENT    0.93 


d.        Regression   Analysis 

To  determine  the  feasibility  of  predicting   mean    sound 
velocity,    a   first-order   (or   linear)    regression   analysis   was   performed   on 
all   data   arrays.      The  data   were  fitted  to  the  equation 

y  =  a   +  bx  (2.20) 

where  y   is   the  estimated   mean   velocity,    a    is   the  coefficient,    b   is   the 
the   intercept,    and    x    is   the  tide   height.      The  coefficient   and    intercept 
estimates   were  determined   to  minimize  the   sum  of  the   residuals 

J^i   =    I^Yi-yj)  (2.21) 

where  Yj  is  the  observed  mean  sound  velocity  and  y-  is  the  estimated 
velocity  [Ref.  13:  p.  533].  Determination  of  the  coefficients  and  the 
intercepts   was   accomplished   by  the   normal    equations 
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an  +  b  2  x:  =    2  Y:  (2.22) 

i=l  i=l     ' 


n  n      o         n 

a    2     X:     +     b    2     X>    =       2     X:Y 

1-1    '         1-1    '         i=l    ' 


(2.23) 


where   n    is   the  number  of  observations  for  a   particular  station   and 
depth    [Ref.    13:    p.    533].      The  sum  of  the  squares  of  the   residuals   is 

SS   =     2    (Yj-Yj)2  (2.24) 

The  standard   error  of  estimate  is 


see  =/(SS  /   v)  (2.25) 

where  v   is   the  degrees  of  freedom  given   &/ 

v  =   p   -    n  (2.26) 

where   n    is   the   number  of   sample  observations   and   p    is   the   number  of 
estimated   variables    [Ref.    14:    p.    131].      For  a   first-order  analysis,    there 
are  two   estimated   variables.       Results   of  a   first-order  analysis   are   listed 
in   Table  VII. 

Confidence   intervals   associated   with    predicted   mean 
velocities   obtained   by   equations   derived   from  the   linear   regression 
analysis   were  calculated   to  determine  the  accuracy  of   predictions. 
Confidence   interval   computation    requires   the   standard   error  of   estimate, 
the   value  of  the   student-t    (t)    distribution   which    is   based   on   the 
desired   percent  of  confidence,    and   the  degrees   of  freedom   for  the 
sample.       For  this   thesis,    90-percent  confidence   intervals   were  computed 
using   the  equation 


CI    =   ±  t(v,0.95)[1    +(1/n)    +  (2.27) 

1 


[(xrx)2  /  2    (x--x)2]]see 
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TABLE 

VII 

4SB--Linear  Regression  —  Depth  =  7  m 

COEFFICIENT  = 

2.8 

Y  INTERCEPT  = 

1468.0 

MEAN 

VELOCITY 

CAST 

RESIDUAL 

1 

-0.3 

2 

0.2 

3 

0.2 

4 

3.2 

5 

2.8 

6 

0.1 

7 

-5.1 

8 

-0.8 

9 

-1.7 

10      a 

-1.2 

11 

-1.0 

12 

-1.4 

13 

0.1 

14 

-0.9 

15 

-4.5  ■ 

16 

-4.0 

17 

1.4 

18 

5.0 

19 

7.2 

20 

5.5 

21 

2.9 

22 

0.0 

23 

-3.2 

24 

-1.2 

25 

0.0 

26 

5.1 

27 

4.8 

28 

0.0 

29 

-1.7 

30 

-0.8 

SS 

244.8 

see 

3.0 

A  RESIDUAL  IS  THE  DIFFERENCE  BETWEEN  THE  OBSERVED 

MEAN  VELOCITY  AND  THE  PREDICTED  MEAN  VELOCITY. 

where  CI    is   the  confidence   interval   for  a   single  estimated   velocity    (y), 
Xj    is   the  tide   height,    x    is   the  mean   tide   height,    and    n    is   the   number  of 
sample  observations    [Ref.    15:    p.    30], 
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Several   plots  of  the  data  arrays   indicated   a   curvilinear 
relationship   between   mean   sound   velocity  and   predicted   tide   height. 
Therefore,    second-   and  third-order  polynomial    regression   analyses   were 
performed  on   all  data  arrays  to  provide  a  better  fit  of  the  data, 
thereby   increasing  the  accuracy  of  the  mean   sound   velocity   prediction. 
The   normal   equation   method   used   in  the  first-order  analysis   was 
extended  for  calculation   of  a   second-order  polynomial.      To  check  the 
results  of  the  programs  for  the   normal   equation   method,    a   subroutine 
using  the  orthogonal   polynomial   method  was  obtained  from  the  NPS 
Computer  System   Library    [Ref.    16:    p.    8].      A  third-order   regresssion 
analysis   was   also  performed   by  the  orthogonal   polynomial    subroutine  to 
investigate  the  possibility  of   increasing   the  accuracy  of  prediction   and 
fit.      Based  on    results  of  the   investigation   and   check    (Chapter   III),    the 
first-order  and   second-order   normal   equation   methods   were  designated 
as   the   primary   methods   of  analysis.      The   normal   equation   method 
appeared   to  be  the  simplest  to  program.      The  second-order  analysis 
uses   the  equation 


y   =  an    +aix    +   a^x 


2  (2.28) 


where  y    represents   the  estimated   mean    sound   velocity   and    x    is   the 
predicted  tide  height    [Ref.    13:    p.    535].      Coefficients    (aQ,    a^  ,    and   a^) 
were  determined   to  minimize  the   sum  of   residuals   by   the  following 
normal   equations 

n  n        0  n 

ann    *   a,    2    x  .  +   a_    2    x.  2  =     2    Y.  (2.29) 

0  li=1     l         2i=i    l  i=1    i 

an  2    x-   a,    2    x.2   *   a„   £    x • 3   =     2    x .  Y  (2.30) 

no  n-3  n/i  n9 

a0.2    x.j2   -   3,2    x.-3   *   a,    2    Xj4     =    .2    x ..  2  Y  (2.31) 

where   n    is   the   number  of  casts   for  a   particular   station    and   depth 
[Ref.    13:    p.    536].      Results   of  a   second-order  analysis   are   listed    in 
Table  VIM. 


31 


TABLE  VIM 

4SB--Polynomial 

Regression  —  Depth  =  7  m 

COEFFICIENTS 

11 

a2 

=  1470.0 

1.6 

0.1 

MEAN 

VELOCITY 

CAST 

RESIDUAL 

1 

-1.7 

2 

0.8 

3 

1.0 

4 

3.2 

5 

2.7 

6 

0.2 

7 

-4*7 

8 

0.6 

9 

-1.0 

10 

-1.7 

11 

-1.4 

12 

-0.9 

13 

-0.7 

14 

-2.4 

15 

-5.0 

16 

-3.2 

17. 

2.0 

18 

4.2 

19 

5.3 

20 

1.4 

21 

-0.6 

22 

-1.7 

23 

-3.0 

24 

-0.4 

25 

0.8 

26 

5.4 

27 

5.0 

28 

0.7 

29 

-0.9 

30 

0.0 

SS 

209.8 

see 

2.8 

A  RESIDUAL  IS  THE 

DIFFERENCE  BETWEEN  THE  OBSERVED 

MEAN  VELOCITY  AND 

THE  PREDICTED  MEAN  VELOCITY. 

Confidence   intervals   for  predicted   velocities   were 
computed    using   Equation   2.27.      The  degrees   of   freedom  for  a 
second-order  analysis   was    less   than   a   first-order  analysis    because  of 
the  addition   of  one  estimated   variable. 
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D.      NOS  BAR  CHECK  DATA 

The  NOS  bar  check  data  consisted  of  46  bar  checks  which  were 
acquired  to  determine  the  corrections  for  soundings  of  survey  H-8422. 
The  format  of  the  data  was   a   listing  of  corrections  for  5-foot   increments 
of  depth.      The  initial   analysis   included  determination   of  the  maximum, 
minimum,    and   mean   correction   for  each   depth    (Table   IX). 


TABLE  IX 

Bar 

Check 

Correction 

Data  For  Survey  H 

-8422 

OBS 

MEAN 
ft 

UPPER 
LIMIT 
ft 

LOWER 
LIMIT 
ft 

RANGE 
ft 

DEPTH 
ft 

88 

-1.0 

-0.6 

-1.4 

0.8 

5 

90 

-0.8 

-0.2 

-1.4 

1.2 

10 

90 

-0.7 

-0.1 

-1.4 

1.3 

15 

90 

-0.7 

.0.0 

-1.3 

1.3 

20 

90 

-0.7 

0.1 

-1.4 

1.5 

25 

87 

-0.7 

0.0 

-1.6 

1.6 

30 

.  82 

-0.6 

0.0 

-1.3 

1.3 

35 

73 

-0.5 

0.0 

-1.2 

1.2 

40 

30 

-0.3 

0.6 

-1.1 

1.7 

45 

12 

-0.4 

0.0 

-1.0 

1.0 

50 

For  each   depth   of  a   bar  check  there  were  two  corrections    listed; 
the  first  correction   was    recorded   as   the   bar  was   being    lowered,    and   the 
second   during    raising.      These  two  corrections   were  averaged   to  obtain 
a   single   value  which   was    used   for  a   correlation   analysis   with    predicted 
tide   height. 
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A  plot  of  the  bar  check  and  tide  height  data   (Figure  2.4)   was 
made  and  visually  examined.      The  plot  indicates  a   poor  correlation 
between   bar  check  corrections  and  predicted  tide  heights.      For 
verification,    a  correlation  analysis    (Table  X)   was   performed   using 
methods  previously  discussed. 


TABLE  X 
Bar  Check  Statistical  Data  for  Si 

MEAN  TIDE  HEIGHT  

irvey 

H-8422 

3  .  3  f  t 

STANDARD  DEVIATION  OF  TIDE  HEIGHT 
MEAN  CORRECTION  

2  .9  ft 

. .  -0. 7  ft 

STANDARD  DEVIATION  OF  CORRECTION  . 

0.3  ft 

CORRELATION  COEFFICIENT  

0.  12 

DEPTH  OF  OBSERVATIONS  

.  .  25.0  ft 

Because  of  the  poor  correlation,    actual   tidal   data   were  obtained   to 
investigate  the  possibility  that  the   results   were   related  to  the   predicted 
tides.      No   improvement   in   the   results   was   accomplished   by   using   actual 
tide  data.      No  further  manipulation   of  the   NOS    bar  check   data   was 
carried  out  because  of  the   number  of   unknowns   associated   with   the 
acquisition    procedures   and   because  the  amount  of   speculation    necessary 
could   not   lead   to  meaningful   analysis. 
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Figure  2.4        Bar  Check   Data   Array--Depth   =   25  ft. 
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III.  DISCUSSION  OF  RESULTS 


A.   CREDDP  DATA 


Data   acquisition   procedures  of  CREDDP  made  an    investigation   of 
mean   sound   velocity  variation   possible,    but  there  are  factors   which 
make  the   study   incomplete.      Because  of  this,    a   complete  model   for  mean 
sound   velocity   prediction   and   correction   based   on   methods   presented    in 
this   thesis  was   not  fully  developed.      However,    the   primary   components 
of  such   a   model   were  devised   and   will   be  explained   in   this   chapter. 
The  factors  which   make  the   investigation    incomplete  and 
recommendations   concerning   verification   of   ideas   and   methods   presented 
in   this  thesis  will    be  addressed   in   Chapter   IV. 

Investigation   of  the  CREDDP   data    revealed  that   sound   velocity    in 
the  Columbia    River  estuary   varies   with   position,    depth,    and   time.      To 
analyze  the  problem  the  study   area   was   divided   into  two  primary 
sections,    the   north   channel   and  the  south   channel    (Table  XI).      The 
'river  mile'   distance   represents  the  shortest  deep-water   path   from   buoy 
"7"   at  the   river  entrance  to  the   station.      The   south   channel,    which    is 
maintained   for  commercial   marine  traffic,    has   a    river  mile  scale  which 
can   be  found  on    NOS   chart   18521,    whereas   north   channel    stations   were 
scaled   using   chart   18521. 

Two  categories  of  mean    sound   velocity   variation   are  considered. 
The  term   "horizontal   variation"   designates   differences    in   mean    sound 
velocity   between    stations;    i.e.,    how  the  velocity  changes   with    location. 
The   horizontal    variation   analysis    is    needed   to  develop   zones   or    regions 
where  velocity   correctors   may   be  grouped.      The  term   "temporal 
variation"   designates   the  change   in   mean    sound   velocity  with   time  at 
each    station.      With   temporal   variation,    mean    sound   velocity   correctors 
for  each    zone  change  with   time. 
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TABLE  XI 

Channel  Data 

north 
Station 

channel                south 
River  Mile         Station 

channel 
River  Mile 

2N 

4NB 

5NB 

2.0             2S 

7.5              4SB 
9.5              5SA 

2.0 

8.0 

15.5 

1 .        Horizontal  Variation 

For  the  analysis  of   horizontal   variation   the  two   requirements 
made  were:    (1)    that  for  each   depth   there  be  but  one   horizontal 
variation   analysis   and;    (2)    the  average  of  all   mean   velocities  obtained   at 
the  depth   being   analyzed   be  designated   as   the  average  velocity.      Use 
of  the  average  velocity  allows   each   station   to   have  a   single  velocity  for 
each   depth.      The  average  velocity   presents   a   problem   if  more 
observations  were  made  on   one  end  of  the  tidal    spectrum  than   the 
other.       For  example,    if  more  casts   were  acquired   during   periods   of 
high   water  than    periods   of   low  water,    the  average  velocity   would   be 
higher  than    if  the  majority  of   casts   were  acquired   during    low   water 
periods.      One  of  the  CREDDP  objectives  was   to  acquire  casts   on   an 
hourly   basis.      This   objective  was   met   in   most  cases,    but  there  were 
some   instances   where  the  time   interval    between    successive  casts   was 
greater  than    1    hour.       For  example,    station   4SB   was   chosen   over   station 
4SA   for  designation   of  the   south   channel    because  there  were   no  casts 
acquired   at   station   4SA   during   the   higher   high   tide.      Another   potential 
problem   would   be   introduced    if  there  were   significant  variations    in   the 
river  flow.      According   to  Jay    [Ref.    5:    p.    114],    river  flow   remained    low 
for  the  acquisition    period.      An   example  of   horizontal    variation   of   the 
average  velocity   is    shown    in    Figure  3.1.      These   plots    indicate  the   rate 
of   horizontal    variation   of  the  average  velocity   is   approximately   the   same 


37 


for  the  north  and  south  channels.      Data  for  the  11-  and   13-meter 
depths  were  limited  and  did   not  cover  the  range  of  tide  well   in   many 
cases. 

Regional   zoning  of  mean   sound  velocity,  correctors   based   solely 
on   horizontal  variation  can   be  accomplished  with   little  difficulty.      Using 
the  depth  with  maximum  horizontal   variation,    zones  can   be  drawn   based 
on  the  NOS  criteria  of  plus  or  minus  4  meters  per  second.      A  major 
problem   still   exists   in   that  the  temporal   variation   of  mean   sound   velocity 
has   not  been   accounted  for  when   only  the   horizontal   variation   of  the 
average  velocity   is   considered. 

2.        Temporal  Variation 

An   analysis  of  temporal   variation   of  the  mean    sound   velocity 
demonstrated   that   proper  correction  *of  soundings   in   the   Columbia    River 
is   a   multi-dimensional   problem.      Sound   velocity   variables   change 
significantly  with   time  at  the  same   location,    and   the  magnitude  of 
variation   can   depend  on   a   number  of  cyclic  and   noncyclic   phenomena. 
Sound   velocity   variables   change  with   the  tide,    and   the  magnitude  of  the 
variation    is   affected   by   monthly  tidal,  events   such   as   spring   and    neap 
tides.      The  temporal   and   horizontal   variations  of  velocity  for  the 
stations   of  the   north   and   south   channels   are  shown    in    Figure  3.1. 

Data  for  stations   2N   and  2S  were  acquired   during   a    neap  tide 
period   when   the  predicted   tide   range  was   minimum,    whereas   data 
acquisition   for  stations  4SA,    4SB,    5SA,and   5SB  was   accomplished 
during   a   spring  tide  period.      Data   for  stations  4NA,    4NB,    5NA,    5NB, 
and   5NC  were  acquired   during   a   period   of  average  tide   range.      As 
illustrated   by    Figure  3.1,    stations   4NB   and  4SB    have  about  the   same 
average  velocity   at  the   standard   depths,    but  with    respect  to  temporal 
variation,    the   range  of  velocity  for  station   4SB   was   significantly   larger 
than   that  found   at  station   4NB.      According  to  Jay    [Ref.    5:    pp 
114,117],    the   neap  tide  period   was   a  one  of  minimum   vertical   mixing   and 
maximum   stratification,    whereas   during   the   spring   tide   period,    maximum 
mixing,    minimum   stratification,    and   maximum   salinity   intrusion   occurred. 
Examination   of  the  casts   for  stations  4NB   and  4SB    (Table  XII) 
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Figure  3.1        Temporal   and   Horizontal   Variations  of  Velocity- 
Depth   =  9  m. 
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confirmed  the  previous  statement.      The  cast  for  station  4NB  acquired  at 
predicted   low  tide  showed  a  significant  difference  in  temperature  when 
compared  to  the  low  tide  cast  of  station  4SB.      The  warmer  temperatures 
of  the  low  tide  cast  caused  the  the  minimum  velocity  to  be  higher  at 
station  4NB  than  at  station  4SB.      Comparison  of  the  high  tide  casts 
indicated  the  surface  layer  salinity  was  6.82  parts   per  thousand   less  at 
station  4NB  than   at  station  4SB.      The   relatively   large  salinity 
difference  for  the   high   tide  casts   caused  the  maximum  mean   velocity  to 
be  lower  at  station   4NB  when   compared  to  the  maximum  mean   velocity   at 
station   4SB.      Therefore,    the  increased   energy  associated   with   a    larger 
tide   range  can    increase  the  mean    sound   velocity   range  for  certain   areas 
in   two  ways.      First,    increased   mixing   associated   with   the   larger  tide 
range  can   moderate  temperature  for   low  tide  cast,    thereby  causing 
minimum  mean   velocities   to   be   lower.      Second,    increased   salinity 
intrusion   associated  with   a   larger  tide   range  can   cause  maximum   mean 
velocities  to  be   higher  when   compared   to  maximum   mean   velocities 
observed   during   an   average  or   neap  tide   range. 

a.  Correlation   of  Tide   Height  and   Sound  Velocity 

The  effort  to  correct  mean   sound   velocity  for  temporal 
variation   began   with   a   correlation   analysis   of  predicted   tide   height  and 
mean   sound   velocity.      The   results    (Table  XIII)    indicate  a   strong 
relationship   between   tide   height  and   mean   sound   velocity  for  most  of 
the  CREDDP  data.      Correlation    in   the   1-meter   layer  was    usually   not  as 
strong   as   the  deeper   layers   for  most  of  the   stations.      Dashed    line 
entries    in   Table  XIII    indicate   insufficient  data   for  a   correlation 
analysis. 

b.  Regression    Analysis 

First-    (linear),    second-,    and  third-order  polynomial 
regression   analyses   were  performed  for  each   standard   depth   of  the 
CREDDP   sampling   stations.      The  objective  of  these  analyses   was  to 
provide  a   numerical   method  of  predicting   mean    sound   velocity   using   the 
predicted   tide   height.      These  methods   would    be  considered   completely 
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TABLE  XII 

The  Effect  of  Tide  Range  on  Mean  So 

und  Velocity 

Range 

STATION   4NB 
TIDE   HEIGHT   =   0.02    ft    ( 

LOW  TIDE) 

DEPTH         TEMP         SAL 
m                °C          g/kg 

VEL 
m/s 

1.0         14.22         6.32 
3.0         14.14         7.48 
5.0         14.05         8.52 
7.0         13.93         9.73 

1472.2 
1472.5 
1473.1 
1473.5 

TIDE    HEIGHT   =    8.3    ft    (HIGH   TIDE) 

DEPTH         TEMP         SAL 
m                °C        g/kg 

VEL 
m/s 

A 

1.0         12.18      23.49 
3.0         11.78      29.26 
5.0         11.97      29.89 
7.0         11.66      30.19 

1484.6 
1487.4 
1489.2 
1489.8 

STATION    4SB 
TIDE    HEIGHT    =    -l.S   ft    ( 

LOW   TIDE) 

DEPTH      TEMP            SAL 
m                °C        g/kg 

VEL 
m/s 

1.0         12.71         8.41 
3.0         12.69         8.92 
5.0         12.66         9.55 
7.0         12.64         9.70 

1468.6 
1468.9 
1469.2 
1469.4 

TIDE    HEIGHT    =    9.2    ft    (HIGH   TIDE) 

DEPTH         TEMP         SAL 
m                °C        g/kg 

VEL 
m/s 

1.3         11.88      30.31 
3.0         11.87      30.69 
5.0         11.84      31.35 
7.0         11.84      31.48 

1491.9 
1492. 1 
1492.5 
1492.7 

successful    if  the  confidence   interval   for   predictions   based   on   all 
possible  tide   heights   was    less   than,    or  equal   to,    plus   or  minus   4  meters 
per   second.      The  employment  of  three  different  orders   of   polynomial 
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TABLE  XIII 

Correlation 

Coefficients 

DEPTH 

STN 

lm 

3m 

5m 

7m 

9m 

11m 

13m 

2N 

0.61 

0.81 

0.85 

0.87 

0.87 

0.87 

0.78 

2S 

0.86 

0.68 

0.79 

0.83 

0.88 

0.89 

0.85 

4NA 

0.83 

0.91 

0.94 

0.97 

0.98 

0.99 

0.99 

4NB 

0.81 

0.90 

0.94 

0.96 

0.96 

0.99 



4SA 

0.71 

0.82 

0.86 

0.89 

0.92 

0.90 

0.93 

4SB 

0.75 

0.83 

0.89 

0.93 

0.96 

0.99 



5NA 

0.86 

0.89 

0.92 

0.91 

0.90 





5NB 

0.80 

0.86 

0.91 

0.94 

0.95 

0.97 

0.96 

5NC 

0.90 

0.91 

0.91 

0.94 

0.96 

0.97 

0.97 

5SA 

0.89 

0.88 

0.88 

0.87 

0.85 

0.86 



5SB 

0.88 

0.88 

0.88 

0.87 

0.86 

- 

regression   analysis   was  for  comparison   of  fit,    as   well   as   comparison   of 
the  predicted   velocities   and   associated   confidence   intervals. 

A   potential   source  of  error  may  exist  when    using 
predictions   calculated   from  tide   height  values   which   are  outside  of  the 
range  of  values   used  for  the   initial    regression   analysis.      For  example, 
if  all   casts   were  acquired   during   a   neap  tide  period,    such   as   was   the 
case  for  stations   2N   and   2S,    velocities   predicted   by  the   regression 
analysis   for   high   or   low  tide   heights  of  a   spring  tide  period   are   subject 
to  error.      For  stations  4SB   and  4SA,    predictions   were  calculated   using 
a  9.0-foot  tide   height;    station   4SB    (Figure  3.2)    was   occupied   during 
higher   high   tide,    whereas    station   4SA    (Figure  3.3)    was    not  occupied 
during  this   period.      The   second-   and  third-order  analysis   predicted 
velocities   which    probably   were   not   physically   possible  at   station   4SA. 
Therefore,    a   basic  objective   in    using   polynomial    regression   models   for 
mean   sound   velocity   predictions   and   corrections   should   be  the 
determination   of  the  maximum  and   minimum  velocities   expected,    and   the 
minimum    range  of  tide  where  these  velocities   could  occur.      The  most 
accurate  method  for  identification   of  these  values   would   be  to  measure 
the  sound   velocity   variables   at  the  same  location   covering   more  phases 
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of  the  lunar  cycle.      The  danger  of  predicting  outside  of  the  observed 
range  is   illustrated   by   Figure  3.3,    whereas    Figure  3.2   shows  that 
predicting  within  the  observed    range  can   be  done  accurately. 

Determination  of  the  most  efficient  polynomial    regression 
analysis  to   use  for  prediction   of  mean   sound   velocity  was   accomplished 
by  examination   of  plots  of  the  data   arrays   and   comparison   of  the 
calculation   of  predicted   velocities   and   confidence  intervals   provided   by 
each   analysis. 

(1)  Data  Array  Plots.        Examination   of  the  data   array 
plots   indicated  that  the  first-order   (linear)    analysis   would   be  sufficient 
for  most  of  the  data   arrays,    although   some  of  the  arrays   did   suggest  a 
curvilinear  trend.      Curvature  could   be  detected   in   the   plots   of   stations 
4NA,    4SB    (Figure  3.2),    4SA    (Figure  3.3),    and  4NB    (Figure  3.4),    with 
maximum  curvature   in   the  plots   of  stations  4SA   and  4SB.      This   may   be 
an    indication   that  the  degree  of  curvature   is    proportional   to  the   range 
of  tide   since  data   for   stations  4SA   and  4SB   were  acquired   during   the 
spring   tide   period,    whereas   data   acquisition   for   stations   4NA   and   4NB 
was   carried   out  during   a   period   of  average  tide   range.      This 
hypothesis   cannot   be  verified   without  further   investigation. 

(2)  Predicted   Velocities .        A   computer   program   was 
written   to   predict  the  mean    velocity  for  five  different   heights   of  tide 
with   the  maximum   and   minimum   heights   being   the   same  as   those   used    in 
the   regression    analysis.      The  other  three   heights   were  equally   spaced 
between   the  maximum   and   minimum.       Rarely   was   there   a   significant 
difference    (plus   or  minus   4  meters   per   second)    between   the   predictions 
calculated   by  the  three  orders   of   polynomial    regression    (Table  XIV). 
Mean    velocities   predicted    using   the  the  third-order   regression    analysis 
were   in    error  for   stations   2S   and   5NB   at  the   13-meter   depth.      The 
cause  of  this   error  could    not   be  explained. 

(3)  Confidence   Intervals .        Confidence   intervals    (Table 
XV)    were  computed   for  the   predicted   velocities    using   t   values   for 
90-percent  confidence.       In    general   terms,    the  confidence   intervals   were 
larger  for  all    stations    in   the   1-   and  3-meter  depths,    with   the   intervals 
becoming    less   with    increasing   depth.      This   corresponds   with   the   lower 
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Figure  3.2       4SB--Prediction   Curves  —  Depth   =   5  m. 
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Figure  3.4       4NB--Prediction   Curves--Depth   =  5  m 
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TABLE  XIV 
Predicted  Velocities—Depth  =  5  m 


STN 
2N 

1.2 

TIDE 
2.6 

HEIGHT  (ft) 

4.0       5.5 

6.9 

ORDER 
1 
2 
3 

1475.5 
1474.8 
1474.1 

1479.0 
1478.9 
1479.3 

1482.5 
1482.7 
1482.5 

•  1486.2 
1486.2 
1485.3 

1489.7 
1489.2 
1488.7 

STN 
2S 

0.3 

TIDE 
1.9 

HEIGHT  (ft) 

3.6       5.2 

6.9 

ORDER 
1 
2 
3 

1476.8 
1479.5 
1479.3 

1480.0 
1479.6 
1479.8 

1483.3 
1481.9 
1481.9 

1486.5 

1486. 1 

■  1486.0 

1489.9 
1492.8 
1493 . 1 

STN 
4NB 

0.0 

TIDE 
2.  1 

HEIGHT  (ft) 

4.1       6.2 

8.3 

ORDER 
1 

2  . 
3 

1472.8 
1474.4 
1474.4 

1476.6 
1475. -9 
1476.0 

1480.3 
1478.8 
1478.8 

1484.2 
1483.5 
1482.9 

1488.0 
1489.6 
1488.5 

STN 
4SB 

-1.5 

TIDE 
1.2 

HEIGHT  (ft) 

3.8       6.5 

9.2 

ORDER 
1 
2 
3 

1465.0 
1469.2 
1469.5 

1471.5 
1471. 1 
1470.9 

1477.8 
1475.9 
1.475.6 

1484.3 
1483.9 
1482.9 

1490.8 
1494.9 
1491.4 

STN 

5NB 

0.2 

TIDE 
2.2 

HEIGHT  (ft) 

4.2       6.2 

8.3 

ORDER 
1 
2 
3 

1468.9 
1470.5 
1471.1 

1472.4 
1472.0 
1471.5 

1476.0 
1475.0 
1475.3 

1479.6 
1479.3 
1480.6 

1483.4 
1485.2 
1485.7 

STN 
5SA 

-1.6 

TIDE 
1.2 

HEIGHT  ( 
4.0 

ft) 

6.7 

9.5 

ORDER 
1 
2 
3 

1460.2 
1461.6 
1461.9 

1463.1 
1462.7 
1462.5 

1465.9 
1465.0 
1464.7 

1468.6 
1468.3 
1467.2 

1471.5 
1473.1 
1468.8 

ORDER 


Degree  of  polynomial  regression 
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correlation  coefficients  found  in  the  shallow  depths.      Larger  confidence 
intervals  were  also  found  for  data  acquired  at  stations  2N,    2S,   4SA, 
and  4SB.      These  stations   had   lower  correlation   coefficients  and   higher 
residuals  when   compared  to  the  other  stations.      For  stations  2N   and   2S, 
this  may  indicate  less   stability  due  to  the  close  proximity  of  the  open 
ocean.      As   stated   before,    data  for  stations  4SA  and  4SB   were  acquired 
during   spring  tide  conditions,    which   tends  to  indicate   less   stability  for 
prediction. 

In   most  cases,    increasing  the  degree  of  the 
regression   analysis   did   not  produce  more  accurate   results   in   terms   of 
confidence  intervals.      The  standard   error  of  estimate  was   sometimes 
lower  with   a   higher  order  analysis,    but  the   reduction    in   the  degrees   of 
freedom  countered  the  effect  of  the   lower   standard   error  of  estimate  for 
computation   of  confidence   intervals. 

Confidence   interval    results   indicate  that  correcting 
soundings  for  mean   velocity  variation    based  on    predictions   provided   by 
the   regression   analysis   cannot  be  considered   a   complete   success. 
Attaining   the  goal   of  predictions   to  within    plus  or  minus  4  meters   per 
second   may   not  be  possible   in    some  cases,    but  the   regression   analysis 
does   show  a   significant   reduction    in   error  when   compared   to  only    using 
average  velocity   and   horizontal   variation   as   the  correction    procedure. 

B.       NOS   BAR   CHECK   DATA 

Analysis  of  the  bar  check   data    revealed   too  many   inconsistencies 
and    required   too  much    speculation   for- comparison   of  the   CREDDP   data 
results.      A   major   problem  of  the  direct  comparison    method   of  the   bar 
check   is   that  examination   of  the   sound   velocity  variables   and   any 
change   in   their  magnitude   is   not  possible.      The  bar  check  data   did 
indicate  the  possibility  of  either  an   extremely   high   sound   velocity 
variation   or  an    instrument  error  which   was   not  constant  with   depth. 
Most   likely,    the   results   show  an    undeterminable  combination   of   sound 
velocity   variation   and   instrument  error.      The  problem  of   using   the   bar 
check   data   for   any   meaningful   analysis    is    illustrated   by   Table  XVI. 
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TABLE  XV 
Confidence  Intervals—Depth  =  5  m 


STN 
2N 

1.2 

TIDE 
2.6 

HEIGHT  (FT) 
4.0 

5.5 

6.9 

ORDER 
1 
2 
3 

4.2 
4.1 
4.3 

3.9 
3.9 

4.0 

3.8 
3.7 
3.9 

3.8 
3.7 
3.9 

4.0 
3.9 

4.  1 

STN 
2S 

0.3 

TIDE 
1.9 

HEIGHT  (FT) 
3.6 

5.2 

6.9 

ORDER 
1 
2 
3 

4.3 
3.8 
3.9 

4.1 
3.6 
3.7 

4.0 
3.5 
3.6 

4.1 
3.6 
3.7 

4.3 
3.7 
3.8 

STN 

4NB 

0.0 

TIDE 
2.  1 

HEIGHT  (FT) 
4.1 

6.2 

8.3 

ORDER 

1 
.   2 

3 

3.1 
2.9 
3.3 

2.9 
2.7 
3.1 

2.8 
2.6 
3.0 

2.8 

'  2.7 

3.0 

3.0 
2.8 
3.2 

STN 
4SB 

-1.5 

TIDE 
1.2 

HEIGHT  (FT) 
3.8 

6.5 

9.2 

ORDER 
1 
2 
3 

5.2 
4.9 

4.9 

5.0 

4.7 
4.7 

4.9 
4.6 
4.6 

4.9 
4.6 
4.6 

5.1 
4.8 
4.8 

STN 
5NB 

0.2 

TIDE 
2.2 

HEIGHT  (FT) 
4.2 

6.2 

8.3 

ORDER 
1 
2 
3 

3.3 
3.0 
2.9 

3.2 
3.0 
2.8 

•  3.2 
2.9 
2.8 

3.2 
2.9 
2.8 

3.3 

3.  1 
2.9 

STN 
5SA 

-1.6 

TIDE 
1.2 

HEIGHT  (FT) 
4.0 

6.7 

9.5 

ORDER 
1 
2 
3 

3.9 
3.6 
4.5 

3.8 
3.5 

4.4 

3.7 
3.5 

4.4 

3.8 
3.5 

4.4 

3.9 
3.6 

4.6 

units   =   ±   m/s 
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TABLE  XVI 

- 

Bar 

Ch 

eck  Correction  Data 

for  Survey  H-8422 

OBS 

MEAN 
ft 

UPPER    LOWER 
LIMIT    LIMIT    RANGE 
ft       ft       ft 

DEPTH 
ft 

DEL  V 
m/s 

88 

-1.0 

-0.6 

-1.4 

0.8 

5 

231.5   • 

90 

-0.8 

-0.2 

-1.4 

1.2 

10 

177.5 

90 

-0.7 

-0.1 

-1.4 

1.3 

15 

124.4 

90 

-0.7 

0.0 

-1.3 

1.3 

20 

93.5 

90 

-0.7 

0.1 

-1.4 

1.5 

25 

88.3 

87 

-0.7 

0.0 

-1.6 

1.6 

30 

78.4 

82 

-0.6 

0.0 

-1.3 

1.3 

35 

53.5 

73 

-0.5 

0.0 

-1.2 

1.2 

40 

44.4 

30 

-0.3 

0.6 

-1.1 

1.7 

45 

55.5 

DEL  V  represents  the  magnitude  of  vel 
necessary  to  compensate  for  the  range 
corrections  if  there  were  no  instrume 

ocity  di 

of  the 
nt  error 

f f erence 
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IV.  CONCLUSIONS  AND  RECOMMENDATIONS 

A.   CONCLUSIONS 

Investigation   of  sound   velocity  variation   in   the  Columbia    River 
estuary   revealed  the  following   significant  findings: 


Sound   velocity  variation    in   certain   areas    required   the  examination 
of   how  veloctiy  changed   with   time  and   location    if   NOS   accuracy 
standards   were  to  be  met.      The  procedure  of  computing  only  one 
corrector  table   based  on   an   average  velocity,    or  correction,    was 
not  acceptable. 

CREDDP  data   analysis   showed   a   strong    relationship   between   tide 
height   and   sound   velocity,    which    stimulated   the  development  of   a 
method  to   account  for  the   variation   of   sound   velocity  with   time. 
Using   a    regression   analysis,    the  accuracy  of  velocity   correctors 
was    improved   significantly,    although   there  were  still   instances 
where   NOS   standards   could   not  be  met. 

NOS   bar  check   data   analysis   showed   a   significant   variation    in 
corrections,    although   the  variation   did    not  correlate  with   tide 
height  as   the  CREDDP  data   did.      Analysis  of  the  variation   found 
at   each   5-foot   increment  of  depth    revealed   inconsistencies   which 
could    not   be  explained. 


B.       RECOMMENDATIONS 

The  following    recommendations   are  made: 


Further   study  of   sound   velocity   variation   must   be  carried   out  for 
verification   of   ideas   and   methods   presented    in   this   thesis. 
Analysis   of  data   acquired   over   longer  time   periods   at  the   same 
locations   will    allow   a   more  complete   study  of   the   variation 
associated   with    different   lunar   phases   and   the   effect  of   changes    in 
river  flow.      NOS   CTD   data   were  acquired    in   the   Columbia    River 
estuary   during   the   same  time   period   as   CREDDP;    these  data   may 
allow   study  ofthe   sound   velocity   variation   over   longer  time 
periods . 

NOS   bar  check   data   analysis   showed   that  the   hydrographer   should 
be  equipped   with   an   efficient  method  of  measuring    sound   velocity 
in   addition   to  determining    instrument  error.      Knowing   the   value  of 
sound   velocity,    detection   of   significant  changes   and   trends    is   more 
efficient  and   accurate.      Bar  checks   should    not   be   used   as   the 

rimary   method  of   correcting   soundings,    but   should   be    reserved 

or  determination   of   instrument  error. 


R 
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Hydrographers  conducting  surveys   in  areas  where  reversing  tides 
could  cause  a  significant  variation   in   sound  velocity  may  benefit  by 
noting  the  height  of  tide  for  each  cast  or  bar  check  acquired  for 
determining  sound  velocity  corrections.      If  significant  variation   is 
detected,    zhe  possibilty  of  predicting   and   correcting  for  temporal 
variation   using  the  tide  height  as  an   independent  variable  should 
be  investigated.      FORTRANT computer  programs  written  for  the 
methods  presented   in  this  thesis  are  listed   in   Appendices  A,    B,    C, 
D,    E,   and   F. 

In   areas  where  temporal   variation  of  sound  velocity  is   significant, 
the  hydrographer  may   benefit  by   using   a   sound   velocimeter  to 
correct  soundings.      The  velocimeter  would    require  only   a   few 
minutes  to  deploy  and    retrieve.      The   hydrographer  could   determine 
sound   velocity   periodically  throughout  the  working   day   without 
significantly   hindering  operations. 
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